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ABSTRACT 

We present the first 3D spectroscopic observations of a nearby Hi detected poststarburst. or 
E+A, galaxy, SDSS J230743.41+152558.4, obtained with the VIMOS IFU spectrograph at ESO 
VLT. Using the NBuRSTS full spectral fitting technique, we derive maps of stellar kinematics, age, 
and metallicity out to 2-3 half-light radii. Our analysis reveals a large-scale rapidly rotating disc 
{vcirc = 300 km s~^) with a positive age gradient (0.6 to 1.5 Gyr), and a very metal-rich central 
region ([Fe/H]=-f0.25 dex). If a merger or interaction is responsible for triggering the starburst, 
the presence of this undisturbed disc suggests a minor merger with a gas-rich satellite as the 
most plausible option, rather than a disruptive major merger. We find spectroscopic evidence 
for the presence of a LINER or AGN. This is an important clue to the feedback mechanism 
that truncated the starburst. The presently observed quiescent phase may well be a temporary 
episode in the galaxy's life. SDSS J230743. 41-1-152558. 4 is gas-rich and may restart forming stars, 
again becoming blue before finally settling at the red sequence. 



Subject headings: galaxies: 
J230743.4H-152558.4) 

1. Introduction 



In the nearby universe, galaxies have a bi- 
modal color distribution, with a blue peak of 
star forming galaxies separated from a red se- 



elliptical and lenticular, cD 



galaxies: individual (SDSS 



quence of quiescent g alaxies b y a sparsely popu- 



lated gap (IStrateva et a l. 2001; lBaldrv et al 



20041: 



Balogh et al.ll2004[ ). The numb er of red galaxie s 
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has roughly doubled since z ~ 1 (jBell et al.ll2004^ . 
Red galaxies, on average, are more luminous than 
blue galaxies. Therefore, a simple cessation of 
star formation in a fraction of the blue galax- 
ies followed by fading and reddening cannot ex- 
plain this build-up of the red sequence. Starbursts 
are often suggested as a way of increasing a blue 
galaxy's stella r mass before letting it fade onto th e 
red sequence (jBell et al.ll2004l: iLabbe et al.ll20"07l) . 
with mergers and interactions as possible starburst 
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triggers (jBekki et al.ll2005l : Idi Matteo et al.ll2008f) . 

Mergers and starbursts have the additional appeal 
that they offer an explanation for the cessation of 
star fo rmation : gas is being consumed by the star- 
burst ( di Matteo et al ] |2007l) an d is being expelled 
by AGN or superno va feedback (jBuvle et al 
Kavirai et al.ll2007[ ). Thus, previously blue galax- 
ies turn red and, through further dry mergers, 
can lead to the formation of the most massive el- 
liptical galaxies observed in the nearby univers e 
(jBell et al.ll2006l: iNaab. Khochfar. Burkertll2006f) . 
This prompted us to investigate the crucial but 
transient poststarburst phase, during which a 
galaxy crosses the color gap, of this particular 
evolutionary pathway. 

Poststarburst galaxies (PSGs, or E-l-A galax- 
ies) have optical spectra with strong Balmer ab- 
sorption lines, revealing young stars, but faint if 
any emission lines, ruling out major star forrnation 




Couch fc SharDleslll987 



Dressier & Gunn|_1198 
Dressier et al. 19991). PSGs offer us a unique win- 
dow on how starbursts affect galaxy evolution. 
The starburst population is fading but still quite 
prominent and we can study both the burst and 
the underlying intermediate- age and/or old stel- 
lar populations. Due to their high surface bright- 
ness, we can study their stellar kinematics out to 
several half-light radii (rg). Clearly, the poststar- 
burst phase is of crucial importance to furthering 
our understanding of the photometric evolution of 
the galaxy population. 

In this Letter^ we present for the first time 
two-dimensional spatially-resolved kinematics and 
stellar population p arameters fo r a gas- rich PSG, 
selected from the iGoto et "aL ( 2003 ) catalog, 
SDSS J230743.41-hl52558.4. With a redshift 
z = 0.0695, luminosity distance = 302 Mpc 
(adopting a spatially flat cosmology with Hq = 
73 km s-i Mpc-\ = 0.27, fl/y = 0.73), dis- 
tance modulus m — M = 37.40 mag, and a spa- 
tial scale 1.27 kpc arcscc"-'^, this is a relatively 
nearby PSG. The few galaxies with known sys- 
temic velocities in a 200 kpc wide search column 
centered on SDSS J230743.4H-152558.4 do not 
coincide in redshift and are spread over a very 
large redshift range (0.036 ^ z < 0.25), mak- 
ing it unlikely that they are physically connected 
with SDSS J230743.41-t-152558.4. We conflrm the 
classification of SDSS J230743.41-hl52558.4 by 
Yagi fc Gotol l|2006h as a field PSG. 



This is the first paper in a series based on our 
VLT/VIMOS observations of 3 PSGs. 

2. Observations and Data Reduction 

The data have been obtained with the inte- 
gral field unit (IFU) mode of the VIMOS spectro- 
graph at ESO/VLT-Mehpal by ESQ staff in mid- 
2006. With the high-resolution "HR-blue" spec- 
tral setup, the 0.67" per lens spatial scale provides 
a field of view (FoV) of 27" x 27" (40x40 lenses). 
This setup offers a spectral resolving power R ^ 
2500 between 4200 and 6300 A. In total, we 
use eleven out of twelve obtained exposures, spa- 
tially dithered in order to work around broken 
IFU fibers, and each with an on-source time of 
2,000 sec. The seeing FWHM varies between 0.7 
and 1.2". Arc line and flat field calibration frames 
have been taken for every Observing Block (OB). 

For the data recuction, we use the generic IFU 
data reduction pipeline implemented in ITT IDL. 
A brief desc ription of the data reductio n steps can 
be found in Chilingarian et al. ( 2007b[ ). Detailed 
information including all VIMOS-specific aspects 
concerning these data will be provided in a forth- 
coming paper. Each of the four VIMOS quad- 
rants, labeled Ql through Q4, of every OB was 
reduced independent from the others until the 
flux calibration step. The relative flux calibra- 
tion of the quadrants is performed using twilight 
sky frames and requires extra care since we use 
spatial dithering. The absolute flux calibration 
was done using observations of spectrophotomet- 
ric standards included in the standard VIMOS cal- 
ibration plan. 

Night sky spectra are reconstructed from the 
lenses outside 5 le of the galaxy to minimize con- 
tamination by the galactic stellar continuum and 
are subtracted from the science frames. Then, 
individual exposures are co-added, applying iter- 
ative sigma-clipping to remove cosmic ray hits, 
spatial and spectral shifts to account for spatial 
dithering, different heliocentric corrections and at- 
mosphere differential refraction. Noise frames are 
computed using photon statistics and are pro- 
cessed through the same reduction steps as the 
science frames, up to the sky subtraction. 

We use twilight and science frames contami- 
nated by Moon (4 out of 12) to assess the data re- 
duction quality and obtain quantitative estimates 
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of the VIMOS spectral line spread function vari- 
ations across the FoV, along the wavelength axis, 
and between different OBs. We fit the solar spec- 
trum, taken from the ELODIE.3.1 stellar library 
(jPrugniel et alj 120071 ). against the observed spec- 
tra in every lens at six wavelength intervals be- 
tween 4000 and 6300A using th e PPXF procedure 
(|Cappellari fc EmsellenJl2004h . From this test, we 
conclude that (1) wavelength calibration quality is 
better than 5 km s'^ between 4000 and 5900A; (2) 
spectral resolution (cj^st) exhibits no significant 
variations across the FoV for individual quadrants, 
neither between observing blocks, although (3) it 
does change along wavelength (60 to 40 km 
from blue to red) and between quadrants; (4) high- 
order moments of the Gauss-Hermite paramet risa- 
tion (/i3 and fe^. lvan der Marel fc Franxlll993l) are 
close to zero except Q3 at A < 4350 A, where is 
modestly negative due to imperfect focusing; (5) 
systematic errors of the wavelength calibration be- 
come very important at A > 5900A, changing sig- 
nificantly between individual OBs so we restrict 
our analysis to shorter wavelengths. 

3. Data Analysis 

We use the NBursts fuU spe ctral fitting 
technique ( Chilingarian et al. 2007a[) with high- 
resolu tion (i? = 10, 000) pegase.hr (|Le Borgne et al, 
20041 ) simple stellar population (SSP) models to 
extract kinematics and stellar population from 
the absorption- line spectra. First, we rebin the 
data to a target signal-to-noise ratio of 15 per 
bin usi ng the Voronoi adaptive t essellation tech- 
nique (|Cappellari fc CopinI 120031 ). We use only 



that part of the VIMOS IFU FoV where the sky- 
subtraction uncertainty is sufficiently small, with 
a size of 18x18 lenses (12" x 12"). The subse- 
quent fitting procedure comprises the following 
steps: (1) a grid of SSP spectra with a fixed set of 
ages (spaced nearly logarithmically from 20 Myr to 
18 Gyr) and metallicities (from —2.0 to -1-0.5 dex) 
is convolved with the wavelength-dependent in- 
strumental response o f VIMOS in every Vorono i 
tessella as explained in Chilingarian et al. ( 2007bl ) ; 
(2) a non-linear least square fitting against an ob- 
served spectrum is done for a template picked 
from the pre-convolved SSP grid using 2D-spline 
interpolation on age (logt), and metallicity (Z), 
broadened with the line-of-sight velocity distri- 
bution (LOSVD) parametrized by v, a, h^, and 



/i4 and multiplied pixel-by-pixel by an n^^ order 
Legendre polynomial (multiplicative continuum), 
resulting in n -I- 7 parameters to be determined 
(we use n ~ 25). 

We exclude narrow 12 A- wide regions around 
the H/3, H7, [Oiii] (A = 4959,5007A), and 
[Ni] (A ^ 5199A) lines from the fit in order 
to be ab le to detect faint e mission lines . As 
shown in IChilingarian et all (|2007bl . I2OO8I ). reli- 
able SSP-equivalent age and metallicity estimates 
are obtained even if the Balmer absorption lines 
are excluded . Usm g the approach described in 
Chilingaria 3 (120091 ) we estimate the excluded re- 
gions to contain only ~ 7% of age-sensitive infor- 
mation. 

4. Results 

We present the maps derived from the full- 
spectrum fitting in Fig [T] Those regions where 
the parameters have large uncertainties (e.g. > 
15 km s~^ for v and a) or are not fit (^13 and /i4 in 
the outer regions) are masked. The isophotes from 
the reconstructed continuum image (datacube col- 
lapsed along the spectral dimension) are shown 
as black solid contours. Line-of-sight velocities 
can be determined with adequate precision out to 
larger radii from the galaxy centre than the veloc- 
ity dispersion and the stellar population parame- 
ters, because they are insensitive to the sky sub- 
traction. 

4.1. Kinematics and dynamics 

SDSS J230743.41+152558.4 exhibits very reg- 
ular, symmetric pure disc rotation with a semi- 
amplitude of 90 km s~^ without any evidence for 
significant isovelae twist. However, there is a sig- 
nificant misalignment of the kinematical and pho- 
tomctrical axes. We mention the evident large- 
scale bar oriented exactly along the galaxy's mi- 
nor axis (roughly North-South) , which is visible on 
SDSS images and also noticeable here in the recon- 
structed image. The rotation curve flattens out at 
~ 4" fro m the centre, corre sponding to 6 kpc or 
- 3.5re (|Buvle et al.ll2006l ) (see Fig. [2|). The ve- 
locity dispersion declines from a central maximum 
of~110kms-i to 40-50 km s-i at 6". The rather 
high v/a in the outskirts of the galaxy also argue 
for the presence of a disc. One has to be care- 
ful with interpreting the high very central values. 
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because a rapidly rotating circumnuclear compo- 
nent smeared by atmospheric seeing may produce 
very similar effects. Weighing all the evidence, 
SDSS J230743.41+152558.4is best classified mor- 
phologically barred lenticular or SBO galaxy. 
This is not surprising since, although the diversity 
is qu ite large, many E +As have SO-like morpholo- 
gies (|Yang et al.ll2008[ ). 

We use a technique si milar to that proposed 
by Ivan Moorsel fc Welld (1985) to determine the 
galaxy's inclination, i, and the position angle of 
the major axis, PA, from the stellar radial ve- 
locity field (the main difference being the non- 
parametric representation of the rotation curve). 
Fitting the velocity field in elliptical annuli in the 
radial range 0.7 < r < 6" gives PA=: 131 deg 
and i = 30 deg. These values agree well with 
the shape and axial ratio of the outer isophote 
displayed in Fig [T] The peak projected velocity 
V ~ 90 km s^^ = W0 sini, after correcting for incli- 
nation effects, yields sa 180 km s^^ for the mean 
tangential velocity. We adopt the value hji = 0.8", 
provided by SDSS, for the exponential scale-length 
of the stellar density profile. Using these num- 
bers, we estimate the circular velocity Ucirc at a 
radius of 4", corrected for asymmetric drift us- 
ing the method outl i ned in paragraph 4.8.2(a) of 
Binnev fc Tremaind ( 2008h . For the outskirts of 
the galaxy, we assume the rotation curve to be 
flat ((70 w (Tij/v^) and the velocity dispersions to 
be independent of radius. We then arrive at the 
following expression for the circular velocity: 



(1) 



Here, ct^ is the tangential velocity dispersion. Un- 
fortunately, along the major axis, cr is a func- 
tion of both (70 and ct^. the vertical component 
of the velocity dispersion tensor. If we assume we 
are dealing with a thin disc then (Tz ^ and 
0-0 ~ a/ sini ~ 80 km s~^. In the case of an 
isotropic velocity dispersion tensor, — az ~ 
cr ^ 40 km s~^. The first case yields fcirc ~ 
300kms~^; the latter leads to Vdic ^ 220kms~^. 
The velocity width (W20) of the Hi 21 cm emission 
of this galaxy is 286±43 km s~^ . This corresponds 
to a circular velocity Wcirc ~ 143/ sin z km s^^ ~ 
286 km s~^ (Buyle et al., in prep.), sitting com- 
fortably between the two extreme values derived 
from stellar kinematics. We will adopt it here as 




our best-guess value. Using these numbers, we es- 
timate the total dynamical mass within a 5 kpc 
radius at - 10" M©. 

With a total B-band absolute magnitude Mb = 
—20.5 mag (Buyle et al., in prep.), logarithm of the 
peak circular velocity log(wcirc) = 2.48, logarithm 
of the central velocity dispersion log((7o) = 2.04, 
SDSS J2307 43.41+ 152558.4 sits exactly on the 
TuUv-Fis her (Il977l) relation for elliptical and SO 
galax i es (iBedregal. Aragon-Salamanca. Merrifield 



I2OO6I : iDe Riicke et al.ll2007t ). It is sl ightly offset 
from the E/SO Faber- Jackson relation ( Faber fc Jackson 
19761 : ICxerhard et al.ll200ll : IPe Riicke et al.ll200,^ . 

being almost ~ 1 mag too bright in the B-band 
for its velocity dispersion. 

I n Fig. [3l we present the Fundament al Plane 
IFF, 



Djorgovski fc Davisll987t ) in K-space ( Bender et al 



1992). On the K2 vs. ki plot (the plane's "face- 
on view" ) , the distinct regions are occupied 
by dwarf dvan Zee et al.ll2004 iGeha et al.l l2003t 



De Riicke et al.l 120051: IChilingarian et all 120081 ) 
and inter mediate luminosity and giant early-type 



galaxies ( Bender et al. 1992) ■ The latter region 



is extended towards the upper-left corner of the 
plot by low-luminosity bul ges and compact ellipti- 
cal g alaxies with M59cO (jChilingarian fc Mamon 
20081 ) being the most extreme case. SDSS 
J230743.41-hl52558.4, indicated by the red filled 
circle, falls in the region occupied by high-surface 
brightness low-luminosity bulges and ellipticals. 
Since the FP is defined for early-type galax- 
ies (E/SO) assumed to be virialiscd systems, we 
should include the rotational energy into the total 
balance at a first approximation as cr^+i'^/S. This 
moves SDSS J230743.41-f 152558.4 to the position 
indicated by the end of the green vector in Fig [31 
toward the locus of bright E/SO galaxies. 

4.2. Stellar populations and nuclear activ- 
ity 

In the stellar population maps, wc clearly see a 
barely resolved core with a shape and size similar 
to that of the central a-bump. The galaxy cen- 
ter reaches a significantly super-solar metallicity, 
up to -1-0.25 dex. The central SSP-equivalent age 
is 650 Myr, increasing to approximately 3 Gyr at 
3" from the galaxy nucleus. The starburst region, 
where the very young po pulation r esides, is about 
2 arcsec (2.5 kpc) across. lYagi fc G oto (20o3) ob- 
served SDSS J230743.41+152558.4 using long-sht 
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spectroscopy and found large H(5 equivalent widths 
out to 2.5 arcsec. The VIMOS data confirm the 
existence of an extended region showing a young 
stellar population, although the positive age gra- 
dient is significant. F or a comparison with the 



SPSS DR6 spectrum (jAdelman-McCarthv et al 



20081 ) of this galaxy, we integrate the light in the 
data-cube inside a 3" circular aperture and apply 
the full-spectrum fitting technique to both spec- 
tra. We find excellent agreement: all the kinemat- 
ical and stellar population parameters, including 
hi, are consistent within the uncertainties. 

We fit the brightest emission line ([Oiii] A = 
5OO7A) with a Gaussian to determine its ampli- 
tude (sec Fig 2^). Unfortunately, no precise mea- 
surements of radial velocities can be made due to 
the very low emission-line fiux. The [Oiii] emis- 
sion is very centrally concentrated, arguing for the 
presence of a spatially unresolved emitting region. 
Careful inspection of the fitting residuals in the 
central region (Fig [JJa) reveals barely-detectable 
emission m the forbidden [Ni] line (A = 5199 A), 
which is a typical signature of a LINER or AGN. 
We therefore connect the central peak of the [Qui] 
distribution to possible nuclear activity. Fitting 
the SDSS DR6 spectrum in the full wavelength 
range of pegase.hr (3900 < A < 6800 A) and 
subtracting the best-fitting template from the 
data reveals strong emission lines in Ha and [Nil] 
(A = 6548,6584A). The high log([Nii6584]/Ha) = 
0.13 and log([Oiii5oo7]/H/3) = 0.46 ratios support 
the LINER interpr etation (jBaldwin et all Il98ll: 
Kewlev et al] l2006h . This would make SDSS 
J230743.41-hl52558.4 the second well-studied 
PSG to date in wh ich low- power nuclear activity 
ha s been detecte d ( Liu et al. 2007.). As proposed 
bv lKavirai et al. (2007), AGN feedback may play 
a crucial role in quenching sta r formation in mas - 
sive PSGs by expcUing the gas. Buvle et al. ( 2008h 
find that most of the Hi gas in the binary PSG 
system EAOIA/B resides outside the stellar bod- 
ies of the galaxies, suggesting a feedback process 
powerful enough to physically displac e large quan- 
tities of gas, such as AGN fe edback ( Silk fc Reed 
19981 : ISchawinski et aLllioOoh . At the same time, 
the faint H/3 and [Olli] (A = 4959A) emission lines 
are detected almost everywhere out to 3" from 
the galaxy centre, suggesting weak ongoing star 
formation over a large part of the galaxy. 



5. Past and future evolution 

The distance to SDSS J230743.41-hl52558.4 
corresponds to a light travel time of ~ 900 Myr. 
We estimate its _B-ban d luminosity evolution us- 
ing the PEGASE.2 (jPioc fc Rocca-Volmerang3 
19971 ) evolutionary synthesis code. If we imag- 
ine SDSS J230743.41+152558.4 evolving passively 
during this period, it will fade by, at maxi- 
mum, 1 mag if all its stars formed 700 Myr 
ago. In the more plausible case of a compos- 
ite stellar population, the effect on the total B 
magnitude will be even smaller, because the old 
population evolves slower than the young one, 
which has a lower mass fraction. This will make 
SDSS J230743.41+152558.4 resemble present- 
day intermediate-luminosity early-type galaxies 
(E/SO) in its position on the Faber-Jackson re - 
lation, the Kormendy relation (|Kormendvlll977^ ■ 
and the Fundamental Plane (see red vector in 
Fig[3]). 4", or 3re, encloses about 75 % of the light, 
resulting in a mass-to-light ratio of M/Lb ~ 4 
in solar units. After 700 Myr of passive evolu- 
tion, this will have increased to M/Lb ^ 10, at 
maximum. This is compatible with pr esent-day 
early-type galaxies (jGerhard et al. l [200ll) . 

SDSS J230743.41-hl52558.4 turns to be a fast 
rotator i f one applies the recen tly proposed classi- 



fication (jEmsellem et al.l 120071 ). However, despite 



its regular morphology, given the presence of a 
dominating large-scale dynamically cold disc, we 
should not classify this object typical early- 
type galaxy. Major mergers and interactions may 
lead to strong, abruptly truncated star formation 
episodes, as sh own statistically using numerical 
simulations by Idi Matteo et al. ( 2007 ). resulting 
in the E-|-A phenomenon. Such violent events 
heat and often even completely destroy the discs, 
leaving little opportunity to explain the observed 
velocity field of SDSS J230743.41-hl52558.4. On 
the other hand, a minor merger of a large disc- 
dominated lenticular or early-type spiral with an 
intermediate-mass gas-rich satellite may be an ac- 
ceptable explanation: young stars will actively 
form in the disc from the accreted intermediate- 
metallicity ISM while the central star formation 
peak will consume the metal-rich gas often present 
in the circumnuclcar regions of early-type galaxies. 



Gi ven the presence of neutral gas (jBuvle et al 
20061 ). the quenching of star formation in SDSS 
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J230743.41+152558.4 could very well be a tran- 
sient phase. For instance, the gaseous compo- 
nent could have been dispersed by supernova or 
AGN feedback and may re-accrete after several 
hundred Myrs, possible restarting the star forma- 
tion. If this object had been observed during a 
star-forming episode, it would have been classified 
as a barred late-type spiral (SBcd). 

Deep Hi radio observations have led to the de- 
tection of 21 cm em ission in about 60 per cent of 
the targeted PSGs (|Buvle et al.ll2006[) . Some of 



the detected PSGs are as gas-rich as local LIRGS 
and mergers, p ossible progenitors o f PSGs in the 



local universe (jKavirai et al.l 120071 ). The pres- 



ence of large quantities of gas indicates that many 
PSGs might actually be observed during an inac- 
tive phase of the star-formation duty cycle. This 
suggests that the evolution from blue to red galaxy 
may encompass several cycles back and forth in 
the colour-magnitude diagram before the galaxy 
finally settles on the red sequence. 
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Fig. 1. — Two-dimensional maps of stellar kine- 
matics, SSP-equivalent age and metallicity for 
SDSS J230743.41+152558.4 derived from the VI- 
MOS IFU data. The black contours indicate 
isophotcs of the reconstructed image, obtained 
by collapsing the data cube along the wavelength 
axis. 




Fig. 2. — The line-of-sight velocity v and 
velocity dispersion a along two cuts through 
SDSS J230743.41+152558.4, one tracing the kine- 
matical minor axis, the other the kinematical ma- 
jor axis. The velocity profile flattens out around 
^ 5" while the velocity dispersion drops to a ~ 
50 km s^^. 
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Fig. 4.— Distribution of [OIII] emission in the central region of SDSS J230743.41+152558.4 (left) and the 
integrated spectrum of the inner 2" (right). 




Fig. 3. — K-space view of the Fundamental Plane. 
The position of SDSS J230743.41+152558.4 is 
shown with a filled red circle. See text for other 
sources of data. The end of the red vector indi- 
cates the position of SDSS J230743.41+152558.4 
after 900 Myr of passive evolution while the green 
vector denotes its shift when taking the rotational 
energy into account. 
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